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Abstract. 

We study the excitonic dynamics of a driven quantum dot under the influence of a 
phonon environment, going beyond the weak exciton-phonon coupling approximation. 
By combining the polaron transform and time-local projection operator techniques we 
develop a master equation that can be valid over a much larger range of exciton-phonon 
coupling strengths and temperatures than the standard weak-coupling approach. For 
the experimentally relevant parameters considered here, we find that the weak-coupling 
and polaron theories give very similar predictions for low temperatures (below 30 K), 
while at higher temperatures we begin to see discrepancies between the two. This 
is due to the fact that, unlike the polaron approach, the weak-coupling theory is 
incapable of capturing multiphonon effects, while it also does not properly account 
for phonon-induced renormalisation of the driving frequency. In particular, we find 
that the weak-coupling theory often overestimates the damping rate when compared 
to that predicted by the polaron theory. Finally, we extend our theory to include non- 
Markovian effects and find that, for the parameters considered here, they have little 
bearing on the excitonic Rabi rotations when plotted as a function of pulse area. 
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1. Introduction 

Semiconductor quantum dots (QDs) provide a promising setting in which to explore the 
interplay of coherent control and decoherence in the solid-state. Spatial confinement of 
charge carriers gives rise to an atomic-like discrete energy level structure within the dot 
region [U EJ E], which allows for the selective probing of particular excitonic (electron- 
hole pair) transitions [H [5] . This has lead to demonstrations of fundamentally quantum 
mechanical effects, such as laser-driven excitonic Rabi rotations [HI [TJ, El El EH EH E21 
[T3l [HI [151 EH] an d two-photon interference in QD emission [TTJ, [HI [191 120]. Moreover, 
optical preparation, control, and readout of a single self-assembled QD spin has been 
achieved [2H [221 [231 [2U [251 I2H1 [27], while various forms of coupling between closely 
spaced dots have been observed and characterised [281 [2S1 EH [311 [32] . 

Such experimental progress clearly demonstrates the feasibility of creating and 
manipulating both excitonic and spin quantum coherence in QD samples. However, 
despite this, QD charge carriers are often still strongly influenced by their surrounding 
solid-state environment. Though the resulting decoherence processes must generally be 
mitigated in order for QDs to be used, for example, in quantum information processing 
devices [331 EH EH EH E7J EH EH HO], they also open up intriguing opportunities for 
exploring system-environment interactions in the solid-state. The combination of strong 
optical-dipole transitions, well developed control techniques, and relatively pronounced 
environmental interactions allows QDs to be used to study important open system effects 
that may be more difficult to observe, for instance, in atomic systems. 

As an example, the damping of excitonic Rabi rotations in single self-assembled 
semiconductor QDs has been demonstrated to be driving- dependent [HI d [HI El EH Ell 
E21E5]. Owing to the large variation in QD type, growth strategy, and experimental set- 
up, a number of possible decoherence channels may be responsible for such behaviour. 
Two prominent mechanisms are off-resonant excitation of the wetting layer [T5 l f¥Tl W2\ . 
and coupling to lattice vibrations (phonons) pEl EH EH Ell ESI EH] ■ In particular, recent 
experiments have provided compelling evidence that interactions with longitudinal 
acoustic (LA) phonons via deformation potential coupling dominates the damping of 
(ground-state) excitonic Rabi rotations in optically driven InGaAs/GaAs QDs [H] [12]. 
Furthermore, while it might be hoped that certain decoherence sources could be 
suppressed by careful selection of samples and experimental techniques, ultimately 
self-assembled QDs are embedded in a host matrix. Interactions with phonons thus 
constitute an intrinsic limitation on the level of coherence seen in their excitonic 
transitions [HI EH EH [46J, H74 EH [49J, EH [51] . As such, a range of theoretical approaches 
have previously been developed to investigate the effects of phonon interactions 
on the coherent manipulation of excitons in QDs. Examples include perturbative 
expansions of the QD-phonon coupling, resulting in master equation descriptions of 
both Markovian [Til \Y2\ H6| [52] and non-Markovian [131 [531 El] nature, correlation 
expansions [HI EH [56] , and non-perturbative, numerically exact techniques which rely 
on calculation of the path integral [45J. 
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The aim of the work presented here is to extend the master equation approach to QD 
exciton-phonon interactions beyond the weak-coupling regime in which it is commonly 
used [57]. The formalism is attractive because in many regimes analytical expressions 
are obtainable, and when they are not, it is not computationally expensive. We present 
theoretical results describing the phonon-induced damping of a resonantly driven QD 
using a polaron transform [37J [53 [59] plus time- local master equation technique [57J 60J . 
Our theory exploits a perturbative expansion in the polaron transformed representation, 
rather than in the system-bath interaction itself. As we shall show, under certain 
conditions this allows us to identify a perturbation term that is small over a much 
larger range of parameters than in the weak-coupling approach [58J. In particular, our 
master equation is able to account for "nonperturbative" effects not captured in a weak- 
coupling treatment, such as multiphonon processes and phonon-induced renormalisation 
of the driving pulse. This is particularly important in exploring the exciton dynamics 
at elevated temperatures (above 30 K for the parameters we consider), where such 
effects may become important. Furthermore, we also extend the master equation to the 
non-Markovian regime. 

We focus particularly on comparing this theory to the weak-coupling Markovian 
technique used in Refs. [HI [12] to provide good fits to both the phonon-induced 
damping and energy shifts of the observed Rabi rotations. We find that for low 
temperatures (below 30 K) the weak-coupling and polaron theories predict essentially 
the same excitonic dynamics, indicating that multphonon and renormalisation processes 
are unimportant. However, as the temperature is increased, we find surprisingly that 
the weak-coupling theory can overestimate the phonon-induced damping rate when 
compared to the polaron approach. This is consistent with the weak-coupling fits to the 
highest temperature plots in Ref. [12]. We also show that, for the laser pulse durations 
used in Refs. [11] and [12] , the inclusion of non-Markovian effects within our formalism 
has an almost negligible influence on the Rabi rotations when plotted as a function of 
pulse area. 

The paper is organised as follows. In Section [2] we introduce a model describing 
the QD system and its interactions with the phonon environment, and also define the 
polaron transformation. Section [3] outlines our master equation derivation together with 
a discussion of the regimes in which we expect it to be valid. In Sections 14. II and 14.21 we 
consider resonant Rabi rotations, driven both by a Guassian laser pulse and constant 
driving, and compare the weak-coupling and polaron theory predictions. In Section 14.31 
we investigate non-Markovian effects, while in Section [5] we give a brief discussion and 
summarise our results. 

2. Model and polaron transformation 

We consider a single QD modelled (as in Refs. [HI [12]) as a two-level system with 
ground-state |0) and single-exciton state \X), separated by an energy oj . The dot is 
driven by a laser of frequency ui, with Rabi frequency ft(t), and is coupled to a phonon 
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bath represented by an infinite collection of harmonic oscillators with frequencies 
and creation (annihilation) operators 6 k (6k). The system-plus-bath Hamiltonian takes 
the form (for h~ = 1) 

H = lo \X)(X\ + n(t) cosooit(\0)(X\ + \X)(0\) 

+ J> k 6 k 6 k + \x)(x\ J2(gA + g*M, (i) 

k k 

where the exciton-phonon couplings are denoted by g k - Moving to a frame rotating at 
the laser frequency ui, and performing a rotating- wave approximation on the driving 
term, we obtain 

Hwa = 6\X)(X\ + ^(\0)(X\ + \X)(0\) 

+5> k 6 k 6 k + \x)(x\ J2(gA + g*M, 

k k 

(2) 

where 5 = uq — u)i is the detuning of the laser from the excitonic transition energy. The 
rotating- wave approximation can be justified here as both the Rabi frequency Q(t) and 
detuning 5 are generally small in comparison to wj. 

To move into the appropriate basis for the subsequent perturbation theory, we now 
apply a unitary polaron transformation to -£/rwa [321 EE] • This transformation displaces 
the bath oscillators when the QD is in its excited state; we shall explore in the following 
section how this can lead in many situations to a smaller perturbation term than a 
weak-coupling treatment of the exciton-phonon interaction in Eq. (jSJ). The transformed 
Hamiltonian is defined by Hp = e s H-^^e' 3 , where 

S=\X)(X\"£(aA-<b k ), (3) 

k 

with «k = fi'k/wk- Hence, we may write 

e ±5 = |0)(0| + |X)(X|n^(±«k), (4) 

k 

where n k .D(±a:k) = e ± ^ k ^ Qkb k _a k 6 k) j s a product of displacement operators £)(±a: k ). 
Defining the Pauli matrices in the (|0), \X)} basis as a x = \X)(0\ + \0)(X\, a y = 
i(\0)(X\ — |X)(0|), and a z = \X)(X\ — |0)(0|, we find that our polaron-transformed 
Hamiltonian reads [581 ES] 

5' ttJt) , n(t) , 

Hp = -a z + —^ L ^x + }^ w k&lA + ~Y~ ( (J xB x + (?yB y ) , (5) 

k 

where the detuning is now 5' = u f — u>i, defined in terms of the bath-shifted QD transition 
energy u' Q = Uq — ^k^kl^kl 2 , and we have ignored an irrelevant term proportional to the 
identity. Bath-induced fluctuations are now described by the Hermitian combinations 

B X = ±(B+ + B--2B), (6) 
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and 

B V = ±(B--B + ), (7) 

where B± = ^Q k D(±«k), and B = (B±) is the expectation value of the bath 
displacement operators. 

Importantly, the driving term in Eq. (j5J) has now been renormalised by a factor equal 
to this expectation value: Q r (t) = Q(t)B. For a phonon bath in thermal equilibrium at 
inverse temperature (3 = 1/ksT, we find 

B = (B±) = exp [ - (1/2) |ct k | 2 coth (/3w k /2)j . (8) 

k 

For models of the type studied here, the system-bath interaction is entirely characterised 
by the spectral density J(co) = J2k \9k\ 2 S(uj — ^k) [HI]- We are specifically interested in 
the coupling of bulk LA-phonons to our QD exciton, shown to dominate the dephasing 
dynamics in Ref. [IT] , and therefore take a spectral density in the continuum limit of 
the form [JTJ [JJ1 H E21 [62] 

J(w) = au 3 e- (u/Wc)2 , (9) 

giving B = exp[— (1/2) J °° du(J(uj)/uj 2 ) coth (/3a>/2)]. The coupling constant a (here 
having units of s 2 ) captures the strength of the exciton-phonon interaction and is 
dependent upon bulk quantities of the QD sample [UJ. The exponential cut-off with 
frequency u c arises from the form-factor of the carrier wavefunctions [HI [62]. For 
excitons in self-assembled QDs it is proportional to the inverse of the carrier localisation 
length, which for simplicity we assume to be the same for both electrons and holes. 

It is important to note that, apart from the rotating- wave approximation on the 
driving, we have made no further approximations in our manipulations leading from 
Eq. (J2]) to Eq. (jSJ) . We have simply put the Hamiltonian into a form that clearly separates 
the effects of the QD-phonon coupling into renormalisation of QD parameters, through 
Q r (t) and oj' , and bath- induced fluctuations, through the last term in Hj>. 

3. Master equation derivation 

Utilising our transformed representation of the QD Hamiltonian, we shall now derive 
a master equation describing the driven QD exciton dynamics under the influence of 
the phonon environment. To proceed, we separate the polaron-transformed Hamiltonian 
such that Hp(t) = H p(t)+Hip(t). Here, H p(t) = Hsp(t)+Hp,p, with bath Hamiltonian 
Hpp = X^k^kfrjc^k, an d time-dependent system part 

H SP (t) = 6 -a z + ^a x , (10) 



while 



Hjp(t) = ^ (a x B x + a y B y ) (11) 
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is the interaction Hamiltonian, to be treated as a perturbation. Moving into the 
interaction picture with respect to H 0P (t) yields an interaction Hamiltonian in the 
(polaron-transformed) interaction picture of the form 

HiAt) = UUt)H lP (t)U 0P (t), (12) 

where U 0P (t) = U SP (t)e' iHBP \ with 



U SP (t) = Texp 



dvHg P (v) 



(13) 



Here, the Schrodinger and interaction pictures have been chosen to coincide at time 
t = 0, while the time-ordering operator T is necessary as, in general, Hs P (t) does 
not commute with itself at two different times. We therefore write the interaction 
Hamiltonian as 

Hj P (t) = ^ (a x (t)B x (t) + & v (t)By(tj) , (14) 

where ai(t) = Ul P (t)aiU SP (t) and B x (t) = e^ BP *^e~ iifBpt , for I — x,y. 

We now follow the standard projection-operator procedure, outlined in Ref. [ST], to 
derive a time-local master equation for the reduced system density operator, psp{t), 
in the polaron frame interaction picture. Considering the QD to be initialised in 
its ground state, with the bath initially in thermal equilibrium, Pb(0) = Pb = 
e -/3 ^ ka ' kfc k bk /tr£(e~ /3 ^ kWkb k fek ) ) we see that the initial system-bath density operator, 
x(0) = |0)(0|pb, is unaffected by transformation into the polaron representation, i.e. 
Xp(0) = e s x(0)e~ s = |0)(0|ps = x(0). Hence, taking a thermal equilibrium state of the 
bath as a reference state and treating Hj P (t) to second order, we find a homogeneous 
equation [53 [60] 

dfcrjt) 
dt 

describing the dynamics of the excitonic system in the polaron frame, under the influence 
of the phonon bath. Substituting in from Eq. (fl4l) . we obtain 
dp SP (t) Q(t) '■' 



/ dstT B ({H lP (t),[H lP (s),p SP (t)p B ]}), (15) 
Jo 



dt 4 



dsn(s) [a x (t), a x (s)p SP (t)}A x (r) 



+[a y (t),a y (s)p SP (t)}A y (r) + R.c.j, (16) 

where H.c. refers to the Hermitian conjugate, and we have made use of the stationarity 
of the bath reference state to write 

(S,(*)A(s))fl = (S,(r)Si(0)>B = A,(r), (17) 
with r = t — s. We note also that owing to the form of B x and B y , correlation 
functions of the type (Bi(t)Bi>(0))b for / ^ V are identically equal to zero. Using 
J(co) = J2k \9k\ 2 $(u — w k) allows us to write the relevant correlation functions in the 
continuum limit as 

A x (r) = ^V« + e -*W-2), (18) 
A,(r) = ^(e^-e-^), (19) 
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where 

r°° / \ 

(20) 



/ duj — r— (cosur coth(/3u;/2) — isinur). 
Jo w v ' 



J drVt{t - t) ([<r x , a x (t - r, t)p SP (t)]A x (T) 



Now, moving back into the Schrodinger picture, and making the change of variables 
s — > t — r, we obtain 

i 

Psp(t) = - ^[S'or z + n r (t)a x , p S p(t)} 
I 

+[a y , a y (t - r, t)p SP (t)}A y (r) + H.c.J , (21) 

where ai(s,t) = Usp{t)Ul P (s)aiUsp(s)Ul P (t). Eq. ( 12T1) is a non-Markovian master 
equation describing the QD exciton dynamics in the polaron frame for a time-dependent 
laser-driving pulse envelope £l(t), and valid to second order in Hip(t). 

3.1. Markov approximation 

While we could directly use Eq. ( I2TT) as a basis for numerical simulation of the exciton 
dynamics (and we shall in fact do so in Section 14. 3f) . a great deal of insight into the 
system behaviour can be gained through the simplifications allowed by the Markov 
approximation. To make a Markov approximation in the present case, we let the 
upper limit of integration in Eq. (l2Tj) go to infinity under the assumption that the 
bath correlation functions Aj(r) decay on a timescale that is short compared to that of 
the system dynamics we would like to capture. Given this, we may also approximate 

E/sp(t-T,t)«exp[i# S p(t)r], (22) 

while replacing Q(t — r) by Q(t) in the integral in Eq. (j2~Ti) . We may then write 

5' 2 cosr]T + Q r (t) 2 S'siwqT 5'n r (t)(l - cost/t) , . 

a x {t-T,t) Rj o~ x -\ a y -\ cr z , (26) 

r] z r] 7] z 



and 



5 sm tit ft r ( t) sin rrr 
aJt — r,t) ~ a x + cosrjTO-y H a z , (24) 



where r) = ^/ 5' 2 + tt r {t) 2 . 

In the following, we shall consider the case of resonant excitation, 5' = 0, which 
simplifies Eqs. ( 123|) and (124")) to a x {t — r, t) = a x and a y (t — r, t) ~ cos(f2 r ()f:)r)(Tj / + 
sm(Q r (t)r)o~ z , respectively. We then arrive at a polaron transformed master equation 

i Q(t) 2 f°° ( 
Psp(t) = ~ -[tt r (t)o- x , p S p(t)) — J dr([a x ,a x p S p(t)}A x (T) 



+ cos(Q r (t)r)[a y , a y p S p(t)]A 



T 



+ sm(n r (t)r)[o-y,a z psp(t)}Ay(r) +H.c), (25) 
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which we shall use to explore the dynamics of a resonantly driven QD beyond the weak 
exciton-phonon coupling regime. For constant driving, Eq. fl25|) is of a familiar Born- 
Markov form in the polaron frame. 

3.2. Regimes of validity 

Having derived Eqs. (!2"T1) and ( J25l) perturbatively in the polaron frame, we should expect 
their validity to be limited in some manner. Recall that the polaron transformation 
displaces the bath oscillators in reaction to a change of state of the QD. Intuitively, we 
would therefore expect the polaron transformed representation to be applicable when 
the bath is able to react on a timescale shorter than, or similar to, that on which the QD 
exciton itself evolves. Since the timescale on which the bath reacts is set approximately 
by the inverse of the cut-off frequency (r B ~ l/w c ), we would therefore expect Eq. (1211) 
to work best in the regime Q/u c < 1. Additionally, the Markov approximation made in 
deriving Eq. ( 1251) limits its validity to timescales greater than tb. 

To put these considerations on a slightly more quantitative footing, we can 
make a rough estimate of the regime of validity of our perturbative expansion by 
considering the magnitude of the perturbative terms in the master equation, namely 
(f2 2 /4)A;(r) [57j [60], for constant driving Q(t) = Q. For example, consider the upper 
bound on the magnitude of A y (r), given by | A y (0) | = (1/2) (1 — B 4 ). Bearing in 
mind that A(r) tends to zero on a timescale of order l/u c , we see that we want 
(fi 2 /4)(|A^(0)|/o; c ) = (Q 2 /8u c )(l - B A ) to be small in the sense that terms higher 
than second order in ifjp may be neglected in the master equation expansion. Since 
(i^ip) = 0, the next term is of fourth order, and its magnitude can be estimated in a 
similar manner by (Q 2 /A) 2 (\A y (0)\ 2 /u>^). Thus, ignoring numerical factors, we find that 
the fourth order term is small in comparison to the second order, provided that the 
condition 



is satisfied [§|. In line with our previous intuition, this condition tells us that in the 
scaling limit {VL/uj c <C 1) we expect our treatment to be valid well beyond a standard 
weak system-bath coupling approach, such that we can explore both the weak 
small a and/or low T) and strong (B <C 1, large a and/or high T) system-bath coupling 
regimes, as well as reliably interpolate between these two extremes (58]. Outside the 
scaling limit, our approach should remain valid provided that the system-bath coupling 
is small enough, or the temperature low enough, such that the inequality of Eq. ( 1261) is 
still satisfied. 

| However, upon transformation back in the original (or "lab") frame changes in the QD state can have 
an influence on the phonon bath. Specifically, the bath state is not stationary in the original frame, 
and the system-bath state is not generally separable. 

§ Considering instead the magnitude of (r) leads to a similar condition, which gives essentially the 
same regime of validity. 




(26) 
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To demonstrate how the polaron approach generally allows a larger regime of 
parameter space to be explored than a standard weak-coupling treatment, we can 
also apply the above reasoning to assess the regime of validity of such a weak- 
coupling master equation. In this case, we have a weak-coupling correlation function 
Aw(t) = J °° dojj(oj)(cos (out) coth ((3u/2) — % sin (ur)) [36], which again falls to zero on 
a timescale of order 1 / u c . Hence, in a similar manner to before, we estimate the second 
order perturbation to be roughly of magnitude |A w (0)|/o; c , while the fourth order is 
then \Aw(0)\ 2 /oJc ■ We then find the condition 



as an estimate of the range of validity of the weak-coupling approach. 

Considering first the zero temperature limit, we find Avy(O) ~ au A , where we use 
the QD spectral density given in Eq. ([9]). Hence, at zero temperature, our condition 
implies that a weak exciton-phonon coupling treatment should be adequate to describe 
the QD excitonic dynamics provided that 



However, as temperature is increased, the magnitude of Aw(0) does too, and we therefore 
expect the weak-coupling treatment to worsen. Approximating coth(/3a;/2) « 2/(/3a;), 
we find 



or au 2 /(f3uj c ) <C 1, which is clearly a harder criterion to fulfill than the zero temperature 
condition. Hence, for a given system-bath coupling strength a and cutoff frequency co c , 
as the temperature of the bath is increased, a weak-coupling treatment of the system- 
bath interaction becomes a worse approximation. 

Though we should be wary of reading too much into numerical values obtained 
from these rough validity conditions, for the system studied in Refs. [HI [12] we can 
take a = 0.027 ps 2 and u c = 2.2 ps _1 extracted through fits to the data, which gives 
au 2 ~ 0.1. Hence, we might expect a weak-coupling treatment to be valid at low 
temperatures for this QD system, as borne out by the excellent agreement between 
experimental data and theory in Refs. jTTl [12]. However, by a temperature of 50 K, 
we find |Avk(0)|/o; 2 ~ 0.4, such that the weak-coupling approximation is now becoming 
dubious. In contrast, for the same parameters, and taking Q = 1 ps -1 , we find that 
the polaron condition [Eq. ([26]) ] gives (Q/uj c ) 2 (1 — B A ) m 0.03 at T = 0, increasing 
up to (Q/u) c ) 2 (l — -B 4 ) ~ 0.15 at T = 50 K. In fact, we shall show below that it is 
around temperatures of 30 K and above that we begin to see significant differences 
between the weak-coupling and polaron treatments of our driven QD, signifying (in this 
case) that the system is beginning to move out of the weak-coupling regime, and both 
driving- renormalisation and multiphonon processes are starting to become important. 




«1, 



(27) 



au 2 < 1. 



(28) 
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4. Resonant excitation dynamics 

We now proceed to explore the excitonic dynamics of our QD system, focusing in 
particular on comparing how the polaron and weak-coupling theories capture the 
interplay between the driving-induced coherent population oscillations and the phonon 
environment as we vary the temperature. Interestingly, we shall see that as the phonon- 
induced damping rate naturally depends upon the renormalised Rabi frequency fl r in 
the polaron theory, but on the original Rabi frequency Q in the weak-coupling theory, 
the weak-coupling approach can actually overestimate the damping rate even in the 
high-temperature regime where mulitphonon effects are important. 

Experimentally, it is generally the excitonic population p X x that is measured, 
for example through photocurrent detection [8j [11] or microcavity-asissted photon 
emission [T3J [H]. Here, we express the solutions to our Markovian master equation 
through the Bloch vector, defined in the polaron frame as ap = (a x p, a y p, a z p) T = 
((a x ) P , (o- y )p, (cr z )p) T , where (a^p = tv s+B (o-iXp(t)), for i = x,y, z. Since a z is invariant 
under the polaron transformation, e s a z e~ s = o~ z , we see that in the original (lab) frame 
ot z = trs+B(o~ z x(t)) = trs+B(o- z xp(t)) = a z p, and the Bloch vector elements along z are 
equivalent in the two representations. Hence, pxx = (1 + ot z )/2 = (1 + a z p)/2, and we 
may work entirely in the polaron frame provided we are only interested in population 
dynamics. 

From Eq. (f25{) we find that the polaron frame Bloch vector evolves according to 

dtp = M(t) ■ a P + b{t), (30) 

where 

/ -(T z -T y ) \ 

M(t) = -T v -a r (t) , (31) 

\ o (ftr(*) + A) -r* J 

and b(t) = (—K x , 0, 0) T . Here, 

r v = ^7.(0), (32) 



2 

n(t) 2 

4 

n(t) 2 



hy(n r (t)) + 7y(-fir(*)) + 27,(0)) , (33) 



where 



and 



A = (S y (flr(t)) - S v {-Qr(t))) > ( 34 ) 

n(t) 2 

K =-y-(7v(nr(*))-7v(-nr(*))), (35) 



7l(w) = 2Re[KP(u)], (36) 
S l (u)=1m[K^\u)l (37) 
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written in terms of the polaron response function 

poo 

kI P \lu) = / dre iujT Ai(r). (3* 



Jo 

Note that as we are solely interested in the exciton population dynamics it suffices to 
consider only the Bloch equations for a y and a z since, in the resonant case, that for ct x 
becomes decoupled. 

The polaron theory developed here will be compared to the Born-Markov weak- 
coupling treatment presented in Refs. [THUS]- With the weak-coupling theory we again 
find that the equation of motion for a x is decoupled, while a y and a z obey 

a y — — T w a y - (Q(t) + \w)®z-, (39) 
a z = Vt{t)a y . (40) 

The damping rate and energy shift can be expressed in our current notation as 

r w = ±( lw {n{t)) + lw (-n(t))), (4i) 

\w = \ (S w (m) - 3 w (-n{t))) , (42) 
respectively, with the weak-coupling correlation function given by 

poo 

A w (t) = / J (u)du! (cos lot coth (/3u/2) — i sinwr). (43) 
Jo 

We can evaluate Y w in closed form, giving 

IV = | J(n(t)) coth(/3fi(t)/2). (44) 

Hence, the weak-coupling rate displays a linear temperature dependence in the high- 
temperature regime |TT] and, as mentioned previously, is dependent upon the original 
Rabi frequency Q(t) as opposed to the bath-renormalised value. Furthermore, we see 
that there is no pure-dephasing contribution to the weak-coupling rate in the Born- 
Markov approximation, in contrast to the terms 7^(0) appearing in the polaron theory 
through and T z . 

We also point out that in the derivation of Eqs. (1391) and (flO]) . no secular 
approximation is made [57]. The corresponding master equation is therefore not of 
Lindblad form and does not guarantee a completely positive, trace preserving map. As 
such, for certain parameters, it is possible that Eqs. (|39|) and (jlQ]) predict unphysical 
behaviour. In fact, we shall see in the following sections that this is indeed the case for 
high enough temperatures. 



4-1. Time- dependent driving 

Having outlined some of the similarities and differences between the weak-coupling and 
polaron transform approaches, we shall now compare their respective predictions in 
the case of resonant driving with a Gaussian pulse envelope. Rather than looking at 
the dynamics in the time domain, we shall instead explore oscillations in the excitonic 
population (Rabi rotations) as a function of varying pulse area, = f_ + ~ n(t)dt, for 
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Pxx 




Figure 1. Excitonic population as a function of driving pulse area (in units of 
7r), for temperatures ranging from 5 K to 75 K, where each curve has been off- 
set by an increasing integer for clarity. Blue solid lines are calculated using the 
polaron approach, while red dashed lines are calculated using weak-coupling theory. 
Parameters: a — 0.027 ps 2 and uj c = 2.2 ps -1 . 



fixed pulse duration, as is common experimentally. We therefore consider a Gaussian 
pulse of fixed width r but varying peak magnitude, centred around t = 0, and described 
by fl(t) = (Q/2Ty/n)exp[— (t/2r) 2 }. Starting at a time —to well before the pulse (i.e. 
t 3> t), we initialise the QD in its ground state: a P = a = (0,0,— 1) T . We then 
numerically solve the Bloch equations (Eq. ([301 in the polaron case, Eqs. ( 1391) and ([401 
in the weak-coupling case) to find the state of the system at any time t satisfying t ^> r, 
such that the pulse has effectively ended. 

Fig. [1] shows the final excitonic population, p X x, calculated from the polaron and 
weak-coupling theories as described above, as a function of total pulse area, (in units of 
7r), for temperatures ranging from T = 5KtoT = 75K (each plot has been offset by an 
increasing integer for clarity). We use experimentally determined values of the exciton- 
phonon coupling strength and cut-off frequency, a = 0.027 ps 2 and u c = 2.2 ps -1 , 
respectively [12], and a Gaussian driving pulse of width r = 4 ps. Note that for the 
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largest pulse areas studied here the ratio fl/u c has a maximum value of ~ 1.2. 

At low to intermediate temperatures (T < 30 K), we see that the weak-coupling 
and polaron theories agree very closely in their predictions for the population dynamics, 
consistent with the excellent agreement found previously between experimental 
observations and the weak-coupling theory in this regime [HI [12] • Importantly, the 
two theories predict almost exactly the same dependence of the Rabi rotation damping 
rate and frequency shift [jj on increasing temperature and pulse area, provided the 
temperature does not increase much above 30 K. As expected, the phonon-induced 
damping is strongly driving-dependent, with oscillations becoming almost totally 
suppressed at high pulse areas for all but the lowest temperatures 0. 

Perhaps the most striking feature apparent from Fig. [TJ however, is that the weak- 
coupling theory tends to overestimate the damping effect of the phonons at higher 
temperatures, when compared to the polaron theory. In fact, as we shall see below 
in the case of constant driving, provided Q/u c ~ 0.7 or smaller, the weak-coupling 
theory predicts a larger damping rate than the polaron theory at all temperatures for 
the realistic parameters studied here. 

At the single-phonon level, this difference can be attributed directly to the 
temperature-dependent suppression of the driving pulse that occurs in the polaron 
transformed Hamiltonian (see Eq. ([5])). Among other things, this has the consequence 
that the rates appearing in the polaron Bloch equations are to be evaluated at the 
(smaller) renormalised pulse strength £l r (t), rather than at the bare pulse strength Q(t) 
as in the weak-coupling theory. The resulting effect can be seen clearly by expanding 
the relevant polaron rates T y and T z (in Eqs. ( |32|) and (|33|) . respectively) up to their 
single-phonon terms. We then find a damping rate of precisely the same form as in the 
weak-coupling theory, 

ri-p h = - J{Q r {t)) coth (f3Q r (t)/2) , (45) 

though evaluated at the renormalised fi r (t), as expected. For low pulse areas, we can 
approximate Fi_ p h w (an/2)Q r (t) 3 coth ((3Q r (t)/2). Hence, for single-phonon processes 
at least, the lessening of the damping rate in the polaron theory is simply due to the 
fact that we are sampling the spectral density at a lower frequency, since Q r (t) < Q(t). 
Any differences would then become more pronounced at higher temperatures, since this 
is when Q r most differs from Q. 

In the full polaron theory, however, the situation is of course much more complicated 
than this simple analysis would suggest. To begin with, we have no particular reason 
to expect the single-phonon rate of Eq. (T4"5"I) to be valid over a larger temperature range 
than the weak-coupling rate of Eq. (jSj), so the sampling of the spectral density at 
different frequencies in the two theories cannot be the whole story. Looking again, for 
example, at the full polaron rate T y (which in fact disappears in the single-phonon 

|| Note that due to phonon-induced frequency shifts the maxima and minima of the curves in Fig. [1] 
are not expected to occur at integer multiples of tt. 

% We do not enter the undamping regime predicted in Ref. [45] for the parameters considered here. 
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approximation), we see from Eqs. (IT51) . (13"2"j) . (I3"6"j) and (TH5T) that it may be written 



Thus, in determining the overall size of the full polaron rates at higher temperatures, 
there additionally exists a competition between the multiphonon effects accounted for 
by the exponentiation of the phonon propagator 0(r), which increases the rate in 
comparison to the single-phonon approximation, and the overall factor proportional 
to Q r (t) 2 , which again acts to decrease it with increasing temperature. 

A further feature to draw out from the comparison presented in Fig. [1] is that while 
the polaron theory predicts physical behaviour at all temperatures considered, for the 
highest temperature (T = 75 K) the weak-coupling theory actually predicts unphysical 
behaviour, since pxx becomes negative for pulse areas ~ tt. This behaviour can be 
related to an overestimate of the phonon-induced frequency shift in the weak-coupling 
analysis at high temperatures, and will again be discussed in more detail below for the 
case of constant driving. 

4-2. Constant driving 

In order to put the arguments outlined in the previous section on a more formal footing, 
it is helpful to consider the dynamics of the QD system for constant driving, in which 
case an analytic form can be given for the population difference, a z = pxx — Poo- We 
construct a second-order differential equation for the time evolution of a z in both the 
polaron and weak-coupling theories. From Eq. (1301) we find for the polaron theory (using 



Y y = n(t) 2 Re / drK x (r) 



n r {tf 



Re / dr(e^ (r) +e" 0(T) -2). (46) 



2 





(47) 



which has solution (for a z (0) 



1) 




(48) 



with damping rate 




(49) 



and oscillation frequency 




(50) 



On the other hand, Eqs. (I39~I) and (140]) give for the weak-coupling theory 




(51) 



which has a solution of exactly the same form, 




(52) 
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Figure 2. (a) Temperature dependence of the weak-coupling rate Tw (red dashed 
curve), polaron rate Tp (blue solid curve), and single-phonon expansion of the polaron 
rate Ti_ p h (green dotted curve), where we have evaluated each rate at ft = 0.5 ps -1 . 
(b) Dependence of the polaron damping rate Tp (blue solid curves) and weak-coupling 
damping rate Tw (red dashed curves) on the driving frequency f2. The two sets of 
curves correspond to temperatures of T — 5 K and T = 75 K, as indicated. Parameters: 
a = 0.027 ps 2 , uj c = 2.2 ps" 1 



though this time with the weak-coupling damping rate IV of Eq. (jHJ), and oscillation 
frequency 

i w = ^4Q(Q + X W )-T 2 W . (53) 

In the constant driving case, we may therefore directly compare the rate Tp and 
frequency £p in the polaron theory to the weak-coupling expressions Y w and £w, 
respectively. 

In Fig. Efa) we plot the damping rates Tp and Tw, along with the single-phonon 
expansion of the polaron rate (T^ph of Eq. (jl5j) ). as a function of temperature for an 
arbitrarily chosen value of the constant driving, Q = 0.5 ps -1 . For these parameters, the 
weak-coupling approximation is indeed shown to predict a larger rate for all values of T. 
Furthermore, there is a significant difference between the full polaron rate (Tp) and its 
single-phonon expansion (T^ph) above temperatures of about 10 — 15 K, indicating that 
multiphonon effects are becoming important. Hence, in this regime, even though the 
weak-coupling rate is still too large, we cannot simply fix it by replacing Q — > Q r in Tw 
(i.e. taking Y w — > r!_ ph ) as this neglects important multiphonon processes. Notice also 
that while the weak-coupling rate varies linearly with temperature above a few Kelvin, 
the single-phonon expansion does not, despite having a very similar form, due to the 
temperature dependence inherent to fl r . 

In Fig. E](b) we show how the polaron and weak-coupling rates vary with the 
strength of the driving frequency Q, for low and high temperatures. As observed 
experimentally [TU [12] , we see a clear and strong dependence on the driving strength 
for both temperature regimes, and in both the weak-coupling and polaron theories. It 
is also interesting to note that around Q/u c ~ 0.7 in the high temperature case, the 
polaron and weak-coupling rates cross, indicating that above this value the hierarchy of 
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Figure 3. Exciton population as a function of pulse area for constant driving: (a) 
using polaron theory; (b) using weak-coupling theory. The different curves in each 
plot correspond to temperatures ranging from 5 K to 75 K in steps of 10 K, with 
lower temperatures coloured blue and higher temperaturea red (the arrows indicate 
increasing temperature). Parameters: pulse duration = 14 ps, a = 0.027 ps 2 and 
lu c = 2.2 ps -1 . 



rates discussed in reference to Fig. |5£a) no longer holds. 

We emphasise again that, as in the case of time-dependent driving, there are 
two important effects present in the polaron theory which are not captured by the 
weak-coupling treatment, and which become increasingly relevant as the temperature 
is increased. Firstly, there are multiphonon contributions, which tend to increase the 
damping rate, as can clearly be seen by comparing the full polaron rate to its single- 
phonon expansion in Fig. E^a). Secondly, the interaction of the QD exciton with the 
phonon bath causes a reduction in the effective driving field. For Q/u c < 1, this tends 
to decrease the damping rate, as can be seen from Fig. [2](b) . 

We are also now in a position to explain the origin of the unphysical behaviour 
predicted by the weak-coupling theory at 75 K (see Fig. [p. In Fig. [3] we again plot the 
excitonic population as a function of pulse area, but this time for a constant driving 
pulse of 14 ps duration, which is roughly equal to the full- width-half-maximum (FWHM) 
of the Gaussian pulse used in Fig. [TJ Notice that for T = 75 K, the excited state again 
takes on unphysical negative values in the weak-coupling theory for pulse areas ~ it 
(right-hand plot (b)). To see how this comes about, we must consider the weak-coupling 
oscillation frequency of Eq. ( 1531) . For > 4f2(f2 + Xw), we find that pxx can take 
on negative values when Xw < — ^, i.e. when the correction to the driving frequency is 
larger than the frequency itself, the weak-coupling theory discussed here breaks down. 
As mentioned previously, this can ultimately be attributed to the fact that no secular 
approximation was made in the derivation of the weak-coupling Bloch equations, which 
therefore have a corresponding master equation which is not of Lindblad form [57] . 

Let us now consider the frequency shift in slightly more detail. In the weak-coupling 
theory we find from Eqs. ( 142]) and (l43|) that 

^^r ^y , (54, 

Jo - w 
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where V indicates that the Cauchy principal value should be taken, and we have made 
use of the identity / °°e* aJs ds = tt5(uj) + V{i/oS). Since we expect the weak-coupling 
theory to break down in the high temperature limit, we can evaluate Xw analytically 
by approximating coth(x) ~ x _1 in the integrand of Eq. (|54|) . In doing so, we find 

X w « Jh£p± (i _ 2(n/ We )F(n/w e )), (55) 

where F(x) = exp[— x 2 ] J* exp[y 2 ]dy is the Dawson integral. The condition X w < — O, 
which determines when we expect unphysical behaviour from the weak-coupling theory, 
then becomes 

< y/nauic, (56) 

where we take the limit Q/u c <C 1. For the parameters of Fig. E^b), we then expect 
to obtain unphysical behaviour when T > 72 K, in good agreement with the actual 
dynamics. We note that while Eq. ( 156]) may give a bound on when the limits of the 
weak-coupling theory are met, its degree of accuracy may become poor well before this 
condition is satisfied (see Eq. ( 129]) and discussion there). 

Turning to the low-temperture regime, where the weak-coupling theory remains 
physical, we have seen previously that by expanding the full polaron damping rate to its 
single-phonon terms we may recover the weak-coupling damping rate, though evaluated 
at a renormalised frequency (compare Eqs. (jUJ) and (H5|) ). In order to complete the 
picture, we shall now show a similar equivalence between the polaron and weak-coupling 
frequency shifts at the single-phonon level. Expanding Eq. ( IM]) to first order in J(co) we 
find the single-phonon approximation to the polaron frequency shift, A — > Ai_ p h, where 

However, this is not quite the whole story since, in the polaron theory, the driving 
frequency is shifted both by A, and also at the Hamiltonian level through O — > O r . 
Ultimately, it is the observables, such as the population difference, that are the physically 
meaningful quantities to consider. Inspection of the frequencies £p in Eq. (150]) and ^w 
in Eq. ( 153]) therefore tells us that we should compare O r (O r + A) in the polaron theory 
to fi(fi + Xw) in the weak-couping theory, as we know that — T y pa Ty/ at the single- 
phonon level. Expanding Q r to first order in J(u), we find Q r (Q r + A) ~ Q(Q + AO) 
where 

AO = O J™ dJM- coth(^/2) ( " a+ ff^~ 1} ) • (58) 

Expanding the remaining factors of B and occurrences of Q r to first order in J(u) we 
find that Eq. (I5B1 reduces to Eq. (15^]) (i.e. AO — > Xw)] m the weak-coupling limit, the 
polaron and weak-coupling theories therefore predict the same correction to the driving 
frequency. 
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t(ps) t(ps) 



Figure 4. Exciton population dynamics in the time domain calculated with (solid 
blue curves) and without (dashed black curves) a Markov approximation. The insets 
show the non-Markovian decay rates (dashed red curves) approaching their constant 
Markovian values (solid orange lines). Parameters: a — 0.027ps 2 , w c = 2.2ps , 
Q = 2ps-\ Temperatures: (a) T = 10 K and (b) T = 50 K. 



4-3. Non-Markovian effects 



Finally, we shall now relax the Markov approximation made in section [37TI to investigate 
non-Markovian effects on the QD exciton dynamics |I3] HH H5j |53j EJ1 [55] [56] , within the 
polaron frame [60]. Referring to Eq. ( f2T|) . and considering the case of constant driving 
for simplicity, we find that our non-Markovian master equation may be written 



Psp(*) 



2 

T 



[5'a z + Q r a x ,p SP (t)) 



dr([a x , a x (t - r, £)p S p(*)]A a 



T 



+ [a y , a y (t - r, t)p S p(t)}A y (r) + H.c). (59) 

Avoiding the Markov approximation thus corresponds to the introduction of time- 
dependent rates and energy shifts in our master equation [57] . 

We determine Bloch equations from Eq. ( 159]) in exactly the same way as the 
Markovian case. Considering resonant excitation, 5' = 0, and inserting a x (t —r,t) = a x 
and o~ y (t—T, t) = cos(Q r T)a y + sin(Q r T)a z (which are exact for constant driving), we find 
an equation of motion for the polaron frame Bloch vector identical to Eqs. (I30 l) -( l35l) but 
with Q(t) — > Q, Q r (t) — > Q r , and all and Si(u) replaced with the time-dependent 
quantities 

-t 



7i (w,t) = 2Re 



'A,(r)dr 



and 



S l (u,t)=lm / e iwr A z (r)dr 



o 



(60) 



(61) 



respectively. For the model we consider here, the difference between the non-Markovian 
and Markovian polaron frame dynamics is entirely captured in Eqs. (1601) and (|6"T1) . The 
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Markov approximation simply corresponds to pushing the upper integration limits to 
infinity. We can therefore make the immediate observation that we should expect the 
Markovian and non-Markovian dynamics to deviate most at short times, since this is 
when Ji(cj, t) differs significantly from 7/(0;, 00) (and similarly for Si(cj, t) and Si(u>, 00)). 
These deviations should be most pronounced when A;(0) is greatest in magnitude, since 
this maximises the difference between the Markovian and non-Markovian rates (and 
energy shifts), and decays on a long timescale (set by l/u c ), as this increases the time 
over which the non-Markovian rates (and energy shifts) reach their Markovian limits. 

To show that this is indeed the case, in Fig. 0] we plot the excitonic population 
of our QD as a function of time (rather than pulse area). For this figure we take the 
relatively large value of Q = 2 ps -1 , so that the excitonic system evolves appreciably 
within the phonon bath correlation time, and in (a) consider a low temperature regime 
of T = 10 K. For these parameters non-Markovian effects are most pronounced at 
short times, as expected, though it is generally fairly difficult to distinguish between 
the two theories, especially beyond t ~ 10 ps. In the inset we plot the non-Markovian 
generalisation of the polaron theory decay rate (see Eq. fj49|) ) 

r P (t) = —{j y (n r ,t) + ly (-ft r ,t)+4j x {o,t)), (62) 

which rapidly approaches its Markovian limit on a timescale ~ 1 ps. 

We can enhance short-time non-Markovian effects by considering higher 
temperatures, as this increases the difference between the Markov and non-Markov 
rates and energy shifts on the bath correlation timescale (though it does not change 
the timescale on which the Markov limit is reached). This is shown in Fig. H^b), where 
we again compare Markovian and non-Markovian dynamics, but now at the higher 
temperature of T = 50 K. Non-Markovian effects are indeed more pronounced at short 
times in this case, and, as shown in the inset, the Markov and non-Markov rates do 
differ more significantly at short times, though the Markov limit is again reached on a 
similar timescale to that at 10 K. Once more, beyond 5 — 10 ps there is very little to 
distinguish the Markovian and non-Markovian dynamics. 

The inclusion of non-Markovian effects within the polaron frame master equation 
can therefore affect the population dynamics at short times (~ 5 ps and below), but 
makes very little difference on longer timescales. When plotting excitonic Rabi rotations 
as a function of pulse area, as in Figs. [T] and [3j it is only the final exciton population 
which is measured. For the parameters of Fig. [3j for example, this corresponds to 
reading out the excited state population after 14 ps. Even for the relatively large Rabi 
frequencies used in Figs. H](a) and (b), we see that non-Markovian effects are almost 
negligible on this timescale. Furthermore, at larger temperatures, for which short-time 
non-Markovian effects seem to be more noticeable, the damping is more pronounced, 
so that the steady state is reached sooner. Hence, since short time behaviour is not 
captured in the pulse area plots of Figs. [T] and [3J neither are non-Markovian effects 
(remember that in Fig. [1] the pulse FWHM is close to 14 ps). In fact, if we plot the 
exciton population as a function of pulse area using our non-Markovian polaron master 
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equation [Eq. (1551) ] for the same parameters as Fig. E](a), we find that it is almost 
indistinguishable from the Markov version on the scale shown in that figure. Hence, 
our theory predicts that (for the parameters considered here) in order for polaron frame 
non-Markovian signatures to be evident in pulse-area plots, FWHM pulse durations on 
the sub 5 ps timescale should be used, much shorter than those in the experiments 
performed in Refs. [TTj IT2] . 

We comment that Markovian and non-Markovian excitonic dynamics could also 
be explored within a time non-local master equation formalism [571 EE]- To do so, 
one would replace the reduced density operator at time t, psp(t) in Eq. (fl5|) . with its 
value instead at time s. This generally results in a master equation of convolution- 
type, which is best solved in Laplace space. The dynamics is then determined by 
the properties of the poles and branch cuts of the Laplace transform of a z (t), which 
will have contributions corresponding to both Markovian and non-Markovian processes. 
Approximating the poles as in Ref. [58J corresponds to including Markovian effects 
only, and should therefore be equivalent to our Markovian treatment. It would be 
interesting to relax this approximation, in particular to include branch cut contributions, 
and compare the resulting dynamics to the time-local approach used in this section, 
since it is expected that both approaches should approximate the true dynamics to the 
same degree of accuracy [57]. We note, however, that a time- dependent driving field is 
likely to cause complications in the time non-local approach, since it would reduce the 
applicability of the convolution theorem. 

5. Discussion and summary 

Inspired by recent experimental observations [TT11T2] . we have investigated the excitonic 
dynamics of a resonantly driven QD under the influence of dephasing due to its 
interactions with an acoustic phonon environment. We have developed a combined 
polaron transform, time-local master equation approach to the problem, which accounts 
for non-perturbative effects such as multiphonon processes and phonon-induced driving 
renormalisation. We have also extended the theory to the non-Markovian regime. We 
have found that for low temperatures (< 30 K), the weak-coupling theory presented 
in Refs. [TT1 [12] is in excellent agreement with the polaron master equation dynamics. 
However, as the temperature is increased, we find that the weak-coupling treatment 
begins to overestimate the damping rate, compared to the polaron theory prediction. In 
fact, it is interesting to note that in Ref. [12] it was reported that a weak-coupling fit to 
the data slightly overestimates the damping for temperatures > 40 K, consistent with our 
findings. For these temperatures, the non-perturbative aspects of the polaron theory 
are becoming important. Renormalisation of the Rabi frequency tends to decrease 
the damping rate, while multiphonon processes act to increase it above the single- 
phonon level (see Eqs. ()44l) . (14"51) and (H9l) for the weak-coupling rate, the single- 
phonon approximation to the polaron rate, and the full polaron rate, respectively). 
Deviations from the weak-coupling theory should be even more pronounced at higher 
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temperatures (above the highest temperature of ~ 50 K explored in Ref. [12]), though 
other decoherence mechanisms could also come into play in this regime. 

We also considered the important role of the energy-shift terms in the weak- 
coupling and polaron theories. These terms, analogous to the Lamb-shift in energy 
levels of atomic physics, are responsible for driving and temperature dependent shifts 
in the exciton population oscillation frequency, as also reported in Ref. [12]. While, 
in general, the energy shifts are necessary for a full description of the dynamics, at 
high temperatures (> 72 K for the parameters studied here) we find that in the weak- 
coupling theory they give rise to unphysical behaviour, and therefore set a bound on the 
applicability of this approach. On the other hand, the polaron theory suffers no such 
limitation in this regime. 

Finally, we explored the role of non-Markovian effects within the polaron frame, 
and found that they are predominantly a short-time phenomenon for our experimentally 
relevant parameters. Hence, they should have little bearing on pulse area plots of Rabi 
rotations if the pulse duration is long on the bath correlation timescale P, El El [H] . Note 
this implies that, under the same excitation conditions, non-Markovian effects are also 
negligible in pulse area plots in the weak-coupling theory at low temperatures, since the 
polaron and weak-coupling approaches agree well in this regime. In order to enhance 
the visibility of non-Markovian effects, shorter duration pulses or longer bath correlation 
times are required. 
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